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Abstract 

The variation of the energy for a gravitational system is directly de- 
fined from the Hamiltonian field equations of General Relativity. When 
the variation of the energy is written in a covariant form it splits into 
two (covariant) contributions: one of them is the Komar energy, while 
the other is the so-called covariant ADM correction term. When specific 
boundary conditions are analyzed one sees that the Komar energy is re- 
lated to the gravitational heat while the ADM correction term plays the 
role of the Helmholtz free energy. These properties allow to establish, 
inside a classical geometric framework, a formal analogy between gravita- 
tion and the laws governing the evolution of a thermodynamical system. 
The analogy applies to stationary spacetimes admitting multiple causal 
horizons as well as to AdS Taub-bolt solutions. 

PACS numbers: 04.20.-q, 04.20. Cv, 04.20. Fy 

1 Introduction 

One of the best-known classical methods to derive the gravitational Hamilto- 
nian from the Hilbert Lagrangian is based on Noether theorem ^61, £3j ; it relies 
on the identification of the Hamiltonian density with the Noether charge density 
associated with a timelike vector field ^. The Hamiltonian is then obtained after 
integration of the Noether current on a (portion) of a Cauchy hypersurface S 
in spacetime. It is also well-known [201 EOl l^^l 1491 1^ that the Hamiltonian 
obtained in this way splits into a volume integral on S (which turns out to 
be vanishing on shell being proportional to the Hamiltonian constraints) and a 
surface integral on 9E which is nothing but the integral of the Noether superpo- 
tential, i.e. the so-called Komar potential |52[ I64| . Therefore the gravitational 
energy, defined as the value on-shell of the Hamiltonian, just reduces to the 
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integral of the Komar form. This simple definition of energy suffers however of 
at least one remarkable drawback. Already when we deal with simple solutions, 
e.g. asymptotically flat stationary spacetimes, it is well known that Noether 
techniques allow one to obtain the expected value for the angular momentum 
but only one-half of the ADM mass. Accordingly, the energy is affected by 
the so-called anomalous factor 1/2 A way to overcome the problem was 

suggested in [301 0H1 (and references therein) and it basically consists in 
correcting the Komar potential by adding to it a new (covariant) background- 
dependent term which allows to cure the anomalous factor. The final expression 
we end up with is nothing but the Noether superpotential associated to the so- 
called first order covariant Lagrangian |EI 12^1 UHl UHl OHI HH ES]> obtained 
from the Hilbert Lagrangian through the addition of a divergence term with an 
arbitrary non-dynamical background dependence. 

The technique of introducing a suitable background dependence into the 
Lagrangian and the superpotential has been shown to work also for spacetime 
solutions other than the asymptotically flat ones and for generic spacetime di- 
mensions (see [IHlliainniEliniEniEniESl) and this naturally leads to the 
interpretation of the corresponding Noether energy as the energy of each solu- 
tion relative to the chosen background.^ 

In this paper, following a thermodynamical-inspired approach |^ EHl 
ISl] . we suggest that the anomalous factor of the Komar potential is, actually, by 
no means "anomalous" . On the contrary what is anomalous is the interpretation 
of Mass/energy we attribute to it. Indeed the Komar potential turns out to be 
the right candidate to describe another kind of gravitational energy, which we 
refer to as the gravitational heat. 

The theory of classical thermodynamics can be of great help in understand- 
ing this issue. A thermodynamic system is endowed with an internal energy 
E. Along a (reversible) thermodynamic process the internal energy can be con- 
verted into mechanical work L and into heat Q according to the first law of 
thermodynamics 6E = 5Q — SL. The amount of energy available to do work 
in a reversible transformation (with initial and final temperature both equal to 
T) is the free energy F, while the difference E ~ F — TS {S = entropy) is that 
part of energy which inevitably transforms into heat; see |27| . 

In the last thirty years a close relationship between the law of thermody- 
namics and black hole mechanics has been established, initially only on a formal 
basis by Bekenstein and Hawking (see [HI El EH|) and subsequently on a more 
physical ground by Hawking's discovery of the black hole radiation; [2^1 OH]- 
Therefore, nowadays we have convincing arguments to state that black holes 
are truly thermodynamic systems (see for a review on the matter). As 
such, we should be able to attribute to them different kinds of energies, such as 
the internal energy, free energy and gravitational heat. According to this per- 
nor course, a number of physically reasonable geometric definitions of gravitational energy 
other than the one based on Noether theorem may be found in literature. Just to mention a 
few of them we recall the Hamiltonian approach and the symplectic methods |^ ^ llbl 1171 1421 
1511 1601 . the Hamilton— Jacobi-based techniques 1101 1131 1141 1151 . and the formulations which 
are directly based on field equations I31II47| . 



2 



spective, the gravitational entropy was calculated in from the gravitational 
heat, suitably defined as the difference between the internal energy (obtained 
via a Brown- York procedure from an action functional with boundary countert- 
erms) and the free energy defined semi-classically as the value on shell of the 
action functional itself. 

In a completely classical (and therefore different) framework we shall give 
here a geometric definition of internal energy and gravitational heat starting 
from a unique generating formula which relates the variation of the Hamiltonian 
with the (pre-)symplectic form of standard General Relativity. Internal energy 
variations SE will be obtained from the variation on-shell of the Hamiltonian. 
When we rewrite 6E in a covariant form it splits into two contributions. One of 
them, namely the variation of the Noether potential, is the only one surviving 
once Neumann boundary conditions on SE are imposed; it then corresponds 
to the Komar energy. The second contribution entering into SE is instead the 
so-called covariant ADM correction term |22 ESI EHl El EH ■ 

When dealing with stationary spacetimes admitting (observer-dependent) 
horizons (hereafter we shall be more precise on the notion of horizon) the Komar 
energy turns out to be the gravitational heat TS (for suitable geometrically 
defined quantities T and S) while the covariant ADM correction term is related 
to the free energy of the system. 

The idea to relate somehow entropy with Noether charges is surely not new; 
see e.g. 0^ (see also where the same idea was worked out within path 
integral methods). We stress however that our definition of entropy differs from 
the original one due to Wald and it is closest in spirit with Mann's prescription 
123 IS] • In this respect we stress that, even if throughout the paper we shall 
deal with the Noether charge, nowhere in the text we shall make in fact use of 
Noether theorem. 

We also stress that the identification of the Komar energy with the gravi- 
tational heat provides just a geometric definition of entropy with no rigorous 
statistical and/or thermodynamical interpretation (even though it turns out to 
coincide with its semi-classically formulation based on path-integral techniques; 
see |^n^^|33). Nevertheless it is supported by a great number of applications. 
Just to mention the most trivial of them let us consider a Schwarzschild black 
hole of mass M. Its temperature T is equal to g^ig- while its entropy S (in 
fundamental units G = h = c = ks = AttM , namely one-quarter of the 
horizon area. Hence the gravitational heat TS turns out to be ^ which is ex- 
actly the Komar superpotential computed on any surface enclosing the horizon 
(and this explains the "anomalous" factor 1/2). This - only apparently - lucky 
coincidence extends to less trivial black hole solutions and, more generally, to 
all stationary spacetimes admitting any number of causal horizons [371 146L I59| 
(including, beyond stationary black hole horizons, also cosmological and Rindler 
accelerated horizons, i.e. horizons which are not associated with black holes). 
Remarkably enough the suggested definition of gravitational heat allows one to 
simply compute the entropy commonly attributed to Taub-Bolt and Taub-nut 
solutions ^1 1^ 1^ 031 ESI Ell whatever the physical meaning is of entropy 
for these solutions. 
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According to our point of view entropy turns out to be a geometric quantity 
the cohomological properties of which can be physically translated into the first 
law of thermodynamics. In agreement with |^ |^ 021 1^^! entropy then 
arises when an obstruction exists to globally foliating spacetime into surfaces 
of constant time. In other words, entropy turns out to be closely related to 
the coordinate singularities of the solution describing spacetime. Being coordi- 
nate singularities a signal that there exist regions of inaccessibility for observers 
which are at rest in the given coordinate frame, entropy eventually turns out 
to be related to the existence of regions hidden to the observers. This perspec- 
tive is close in spirit with Padmanabham's recent and interesting work relating 
entropy with unobserved degrees of freedom; see j^Hl • In agreement with Shan- 
non's original formulation of the theory of information, entropy encodes thence 
the information content hidden beyond inaccessible regions. Not surprisingly 
our definition of entropy agrees with the definition given in |59j for static space- 
times, thus endowing Padmanabhan's physically-motivated construction with a 
mathematical ground and a much broader domain of applicability. 

The paper is organized as follows. Starting from Hamilton's field equations, 
in section 121 we derive the relevant expression for the variation SE of the energy. 
We rewrite 6E into an explicit covariant form and, afterwards, we focus our 
attention on boundary conditions. Depending on the choice made between 
Dirichlet or Neumann boundary conditions we obtain different definitions of 
quasilocal energies, i.e. Dirichlet energy or Neumann energy. The physical 
interpretation of such boundary conditions is carried out in sectional It is here 
recognized that, in spacetimes admitting causal horizons, the Komar energy 
is related to the gravitational heat, while the ADM correction term plays the 
role of the free energy. Moreover, the cohomological properties the variational 
formula for SE is endowed with allow one to relate the entropy to homological 
obstructions in spacetime. We generalize this idea in section^ We suggest, in 
complete agreement with [SJ , that entropy arises whenever an obstruction exists 
in foliating spacetime into space -I- time, thus corroborating previous results 
with an alternative mathematical viewpoint. Moreover, the interplay between 
boundary conditions and the choice of a reference background is discussed. As 
a good example, in section we analyse the AdS Taub-bolt solution and we 
reproduce within our formalism the same numerical values found elsewhere by 
means of different techniques. 

The appendix deals with spherically symmetric solutions with a cosmo- 
logical constant (Schwarzschild-deSitter) or with generic external matter fields. 
This appendix is introduced in order to test on a simple model the main formu- 
lae introduced throughout the paper, thus endowing them with a direct physical 
interpretation. Finally, in appendix ^ we carry out the analogy between the 
classical geometric framework developed in this paper and the semi-classical 
statistical approach based on path integral techniques. 
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2 The Generating Formula and Boundary Con- 
ditions 



Let us consider, in a Lorentzian spacetime M of dimension 4, a region D that is 
difFeomorphic to the product T, x IR where S is a 3~dimensional manifold with 
a sufficiently regular boundary B = dT,. ^ We denote this diffeomorphism by 

: S X J? — >D (1) 

For any i e a hypersurface C I? is induced by ip according to the rule 
Et = ^/'(S X {t}) = V't(S) and we require it to be a (portion) of a spacelike 
Cauchy hypersurface. The set of all St, for t in M, defines a foliation of D la- 
belled by the time parameter t. Moreover, each Et intersects the boundary dD 
in a 2-dimensional surface Bt which is diffeomorphic to B, for all t in M. The 
diffeomorphism is established by the map ipt ■ B — * Bt. Hence the boundary 
B = dD is a timelike hypersurface globally diffeomorphic to the product mani- 
fold B X R and, physically speaking, it describes the histories of the observers 
located on Bf . The time evolution field ^ in _D is defined through the (local) 
rule f Vpt — 1 and, on the boundary B, it is tangent to the boundary itself. 
We denote by the future directed unit normal to Ej and we denote by 
the outward pointing unit normal of Bt in E^. For simplicity we assume the 
foliation to be orthogonal so that the vector is also the unit normal of B in 
spacetime M. Accordingly, everywhere on B it holds true that u^n^jg = 0. 
The evolution vector field can be decomposed as 

^^ = iVw^ + 7V^ (2) 

where N is the lapse and the shift vector iV^ is tangent to the hypersurfaces 
Et, i.e. Nf^ Uf_, — and N^^ n^jg = 0. The metrics induced on Et, B and Bt by 
the metric are given, respectively, by: 

= + Ui, (3) 

7^1^ = 9fj^i' - n-f^ iiu (4) 
o-pi. = h^_,u - ni_, ny = 7^^. -|- u^, (5) 

The metrics Q, Q and Q with an index raised through the contravariant 
metric g^^" define the projection operators in the corresponding surfaces. We 
also denote by 

/i^V„u. (6) 
7^V„n, (7) 
a^D„n, (8) 

the extrinsic curvatures, respectively, of Ej in Af, of B in M and of Bt in E(. The 
symbol D denotes the (metric) covariant derivative on E^ compatible with /i^j,. 

■^Even though the formalism can be generalized to any spacetime dimension > 2, for the 
sake of simplicity we assume dim M = 4. 
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We shall denote hy K = K^^h'''', 9 = e^^7^'' and K = /C^^cr^^ the traces of 
the appropriate extrinsic curvatures. The momentum P^^'^ of the hypersurface 
Et is defined as: 

while for the m,omentum li^'^ of the hypersurface B we set: 

n^" = (67^^ - e^") (10) 

Let us consider the projection of Einstein's equations (in vacuum or with a 
cosmological constant) into their normal and tangential components relative to 
the leaves St of the foliation: 

H = Q (11) 

Ha = (12) 

£$,h^v = [hf_,u] (13) 
_£^p/x-^ = [p/^-^] (14) 

where Ti,Tia are the usual Hamiltonian constraints [/i/ijy] — 2^(2P^i^ — 
h^vP) + "^-Di^^Ni,) and [P^"] shortly denote the right hand side of the Hamil- 
ton equations, the detailed expression of which is of no primary interest here 
and can be found in jHZ]. Notice that is nothing but the definition of 
the momentum P^'^ coniugated to the spatial metric hfj_^, i.e. the Legendre 
transformation, while equations H12|l and 1)14(1 correspond, respectively, to 
Einstein's equations G^'^u^u^, G^^'^ha^Ui, and G°'^h'^h'^ written in a first-order 
formalism in terms of the conjugated fields h^^i, and P'^'^. 

Let us now denote hy X = (5g^^ a vertical vector field in the configu- 
ration bundle of the Lorentzian metrics (in the sequel we shall shortly denote 
with Sy the variation induced on a generic metric-dependent field y by the 
variation X of the metric itself). If we multiply the equations pi|) - (|14|l with, 
respectively, SN, SN", 5P^^ and Sh^^, we sum up all the obtained expressions 
and we suitably rearrange the terms, we finally obtain that the Hamiltonian 
generating equations I(1HI - ((14|I is constrained by the variational formula: 



SxHi^,j:t)^^i^t,X,£^g) (15) 



where: 



SxH{^,j:t)= f {HSN + n^,SN^ + [h^,]SP^'' -[P'''']Sh^,}d''x (16) 



I St 

denotes the variation of the Hamiltonian along X, while 



(17) 
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is the (pre-) symplectic form; see After a lenghty calculation (see, e.g. 

^lE]), expression H16() can be conveniently rewritten as: 

5xHiC,^t) ^ Sx [ {Nn + N^Hajd^x (18) 
+ j^^ d^x |iV^(V^e) - 7V"<5(V^j„) + ^^s'^Ha^p^ 

where 

(20) 

(a^ = u^V^w^) (21) 

are, respectively, the surface quasilocal energy^ the surface momentum and the 
surface stress tensor which describe the stress energy-momentum content of the 
gravitational field inside Bt (see ^3 El E] and references quoted therein) . 
We point out that for non-orthogonal foliations of spacetime, i. e. w^n^jg ^ 
0, an extra term has to be added in IjlSfl which corresponds to a symplectic 
boundary structure; see [TUIITBI 1551 I5T] . 

We define the (variation of the) energy contained in a region St enclosed by 
a surface Bt as the value on-shell of the (variation of the) Hamiltonian. From 
H18|l we obtain: 



5xE{^,Bt)^ 




where now the vector field X — Sg^^-^^ is meant to be a solution of the 
linearized field equations or, equivalently, it can be viewed as a vector tangent 
to the space of solutions (which in turn can be identified with the phase space 
once the suitable gauge reductions have been performed; see JS|)- Moreover 
from H15|) we have the on-shell relation: 

5xE{S„d^t)^ I Lu{X,£^g) (23) 

Hence, for stationary spacetimes which satisfy vacuum Einstein's field equations 
(possibly with a cosmological constant), the 2-form that defines the variation 
of energy is a closed form. Indeed, the right hand side of H23|l is vanishing when 
^ is a Killing vector: 

SxE{C,dEt) = if £e.g = (24) 

We stress that, once we fix Bt in the expression (|22ll . there exist many 
"energies" , depending both on the symmetry vector field ^ and on the variational 
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vector field X; see, e. g. P El |7| El ESI El ■ Tlie former vector determines, via 
its flow parameter, how the observers located on Bt evolve as the "time" flows. 
The vector field X fixes instead the boundary conditions on the dynamical fields 
(e.g. Dirichlet or Neumann boundary conditions). Each choice of the pair (^, X) 
gives rise to a particular realization of a physical system endowed thence with 
its own "energy content" . 

Remark 2.1 We stress that, despite the fact that the Hamiltonian formulation 
breaks down the covariance of the theory through the spacetime foliation into 
space+time, covariance can be nevertheless restored at the end. Indeed formula 
H22f) can be conveniently rewritten in an explicit covariant form as follows: 

5xE{^,Bt)^ j ]-U''P{^,X)dsc.f3 {ds^f3^ip\io,\d^x) (25) 

J Bt 



where 



+ ^5''''H^.rl (26) 



with u^jj = rj^^ — (Jl'^r^j^. The first term in the right hand side of (|26|l is nothing 
but the variation of the Noether superpotential ?7Kom(0i namely the Komar 
potential, while the second term is the so-called covariant ADM correction term 

C/cadm(C,^) (see Q2.2^MW)- 

Written in terms of variables adapted to the foliation they read, respectively, 
as follows HSlEni^: 

5x1 UKon.{0 = ^ [ d^xS[2V^u^e';^^ (27) 

JBt ^'^ JBt 



and 

t^CADM(^, X) 

Bt 



- I d^x^^,5W 

JBt 

!'« Jst 



(28) 



Taking into account the useful formulae 1141 : 

©Miy = + {n°'aa)u^u^ + 2a'^^u^)Kai3n'^ (29) 
(57^, = ~(2/N)u^uJN - (2/iV)(7„(^u,)(57V" + 'j'^y^-^Sa^p (30) 

it easily follows that the sum of (|27|) and H28|) corresponds to (|22|l . 
We just stress that it is indeed the sphtting 

5xE{i, Bt) = 5x C/Kom(C)| + C/cadm(^, X) (31) 

•^We point out that, if we had not assumed the foliation to be orthogonal, there would 
appear additional terms in the right hand side of formulae I2YI and 1281 which however vanish 
when we assume J to be a Killing vector. With the assumption of orthogonal boundaries these 
extra terms are nevertheless vanishing without any Killing requirement; see 1381145) . 
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which will come naturally into play when attempting to describe the gravita- 
tional system in term of thermodynamical variables. Indeed the two contribu- 
tions in (|31|l will be identified, respectively, with the S{TS) and SF contributions 
to the internal energy E (where F denotes the Helmholtz potential) . 

We also stress that was obtained in ^H^HIIS^^ in a purely Lagrangian 
framework as the variation of the (covariantly conserved) charge associated to 
the generator of symmetries ^. From the explicit covariance of H31|l it turns 
then out that the variation of energy depends on the generator of symmetries 
^, namely it depends on the observer, but it does not depend on the coordinate 
system. ■ 

We point out that the expression H22(l can sometimes be integrated in specific 
applications, providing us with an explicit value of E; see e.g. section El and 
Appendix Nevertheless there is no guarantee in general that there really 
exists a function E such that its variation fulfills H22|l . even if the infinitesimal 
quantity SxE is always well defined in our hypotheses. 

The problem we address now is to implement boundary conditions which 
allow formal integration of expression 1)22(1 . 

Dirichlet energy. Wc define the Dirichlet quasilocal energy £{£^, Bt) contained 
in the surface Bt to be the value obtained from l|22|l by imposing the Dirichlet 
boundary conditions Sjfj_^\^ ~ 0. Since (57;^i/|g = implies (5A^|g = 0, (JA^^jg = 
and J(T^i/|g = 0, formula (|22|l can be explicitly integrated: 

£{^,Bt)= f d^xV^{Ne~N''jc.)+£o (32) 

JBt 

We stress that we are integrating in the space of solutions along a curve of solu- 
tions (i.e. the integral curve of X through the point g) satisfying the boundary 
Dirichlet conditions. Therefore the constant of integration £q can be considered 
as the quasilocal energy associated to a point go in the curve and can be fixed 
as the reference point (or zero level) for the energy namely 

£{tBt)= ( d^x^{N{e~eo)-N"{3c.-joc.)} (33) 

JBt 

where the subscript refers to the background solution g^. Comparing ((83|l and 
(|22(l we obtain the relevant formula 

5xE{^, Bt) = 6x£{^, Bt) + f n^"^ 57^, d^x (34) 

JBt 

which obviously reduces to S^jj^E ~ (5[_d]£ when we consider variations generated 
by a vector field X satisfying the Dirichlet boundary conditions <^[_D]7piy = 0. 

We remark that formula (|32|) is nothing but the value on-shell of the Hamil- 
tonian ensuing from the canonical reduction of the Trace-K Lagrangian |15[ I25L 

EH EH. 
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When the observers located on Bt evolve with velocity = u'^(i.e. N — 1, 
iV^ = 0), the Dirichlet energy 1)33(1 reduces to the Brown- York quasilocal energy 
£by{£,, Bt) = d?x^/a{e — cq); see ^i]. The relationships between the energy 
H33|l and the Brown- York energy have been analysed in ^1 ^1 ^] (see also 
|25|). We also point out that there is not a commonly accepted and preferred 
choice for the (variation of) energy of a gravitational system among 5x£{£,, Bt), 
Sx^Bvii, Bt j) or dxE{^,Bt). We nevertheless stress that both the Dirichlet 
expression Sx£{^,Bt) and the Brown- York definition Sx£by{£,, Bt) are a par- 
ticular case of the more general prescription (|22() . In the sequel we shall make 
use of this latter definition. Indeed, we remark that in all the specific solutions 
analysed so far (including Kerr-Newman, BTZ, Taub-bolt, isolated horizons 
solutions) expression 122)) always gives rise to the expected numerical value; see 
PI23liniE3| and the appendix |^ 

Neumann energy. In deriving expression ()32)l we have selected the compo- 
nents of the boundary metric as the symplectic control parameters of the 
gravitational system |17l 151) . while the response variables are the components 
of the boundary momentum H'^'', namely the quasilocal surface energy y/ae, the 
quasilocal surface momentum y^ja and the surface pressure y^s"^ . 

It was argued in (121 that boundary conditions in General Relativity 
exactly correspond to boundary conditions of a thermodynamical ensemble and 
that dynamical fields which are conjugated to each other in a symplectic sense, 
can be also considered as being thermodynamically conjugate. Accordingly, 
different ensembles can be realized by exchanging the fields that are kept fixed 
as boundary data with their respective boundary conjugated fields. Under this 
point of view the variational vector field X turns out to be a variation into the 
ensemble realized with specified boundary conditions. 

In obtaining formula ()32)l we kept fixed the boundary 3-metric and we ob- 
tained the Dirichlet energy. The boundary fields which are symplectically con- 
jugated to the boundary metric are the boundary momenta 11'^'^ defined in ()10)) . 
Hence we formally describe the gravitational system in another ensemble by 
taking the momenta as boundary fixed data. We shall now check if the new - 
Neumann - boundary conditions lead to a mathematically well-defined notion 
of energy, i.e. we will check whether formula 1)22)1 turns out to be integrable 
under the conditions (511^" |g = 0. That this is indeed the case easily follows 
by observing that the covariant ADM potential ()28)l is vanishing if the much 
weaker condition 

-f^JW\^ = (35) 

is imposed and, consequently, the contribution to the variation of energy comes 
entirely from the Komar potential, i.e. 

<5[Ar]i?(e,i?t) = Vl / C^Kom(0 (36) 
JBt 

where now S[n] refers to variations generated by a vector field X satisfying the 
(weak) Neumann boundary condition ()35)l . Therefore we end up with a new 
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energy, say the Neumann quasilocal energy Q, defined as 

QitBt) = I UKon. + Qo 

^ - I ^(fx{Nn^a^~N''K^pn!^]+Qo (37) 

J Bt 

where in the second equahty we made use of (P7|l and (j^ . The constant 
of integration Qq again refers to the energy of a (background) solution lying 
in the integral curve of X through the point g. 
Comparing now H37|l with (|22|l we have: 

5xE{i, Bt) = SxQi^, Bt) - [ 7m- '^n'"' d'^x (38) 

JBt 

Finally, from (|34(l and (|38(l we infer that 

6xS{t Bt) = Sx |g(e, Bt) - n^"^ d^xj (39) 

i.e., the Dirichlet quasilocal energy £ and the Neumann quasilocal energy Q are 
related by a boundary Legendre transformation exchanging boundary metric 
with boundary momenta. 

We shall now analyse to what extent we can refer to the energy (|37|l as the 
physical gravitational heat of a gravitational system. 

3 Entropy of Causal Horizons 

In this section we shall deal with 4-dimensional stationary spacetimes admitting 
one or more causal horizons; see |37l 1461 1^ . A causal horizon Ti is defined as 
the boundary of the past of a (future inextensible) timclikc curve which repre- 
sents the observer's worldline. This observer dependent definition encompasses 
as a special case all (observer independent) standard black hole horizons but it 
also includes observer dependent horizons such as the cosmological horizon as 
well as the Rindler horizon (i.e. the horizon tested by an accelerated observer in 
Minkowski spacetime). In a recent paper ^S] the black hole laws of thermody- 
namics were extended to the more general setting of causal horizons by relying 
on the notion of local horizon entropy. We claim that the same local notion of 
entropy can be formulated in terms of the Komar energy (|37|l . extending in this 
way Wald's original definition IHJ to a broader context. 

In the previous section we gave the definition (|22|) for the (variation of the) 
internal energy contained in a spacelike region Ef bounded by a 2-dimensional 
surface Bt- Let us now assume that the spacetime is stationary and that the 
region contains the cross section H of a, causal horizon Ti. Generalizing 
Wald's original definition for stationary black holes, the horizon entropy S can 
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be heuristically defined for any H as that quantity satisfying the first law of 
thermodynamics: 

5xEitH)^TSxS (40) 

If the temperature T can be provided a priori by physical arguments, thence 
the function S can be truly identified a posteriori with the physical entropy. 
Namely, expression (|40|l can be considered as an algorithm to define a quantity 
S which, by its own definition, is related to the variation of energy in order that 
the first law holds true. 

Let us therefore recall what a physically-motivated definition of tempera- 
ture could be. In |37| it was shown that causal horizons in stationary spacetimes 
which satisfy Einstein's equations in vacuum or with a cosmological constant 
(and possibly with an electromagnetic field) are necessarily stationary axisym- 
metric Killing horizons. Let us then identify the vector field (|2l with the vector 
field which coincides on the horizon with the null geodesic generator (if multiple 
horizons Tit are present, each one is endowed with its own vector and two 
different vectors differ on any horizon for an axial Killing vector, the orbits of 
which are closed curves of length 27r; see |23E1)- In analogy with the black 
hole case, the temperature T of a causal horizon Ti. is identified with the surface 
gravity, i.e.: 



27r 



477 



Nnf^afj, 



2tt 



(41) 

H 



H 

which is constant on Ti. |46| and agrees with the Hawking temperature for black 
hole horizons [35] or cosmological horizons ^7\. For spacetimes with a single 
horizon, the (global) temperature can be equivalently defined as the inverse 
of the period (3 of the Euclidean time. Namely, when we approach closely to 
the horizon along a radial direction (let us suppose that the metric becomes 
spherically symmetric as it approaches the horizon^) the complexified metric 
can be written in its Rindler form: 

^2 C^^V ,2 



g^R' j dr' + dR^ + r^{R)dn (42) 

where t = it is the complexified time, dQ is the metric of the unit sphere and 
/3 — 1 /T is given by 1)41(1 . Moreover 

r{R) = ^R^ + r+ (43) 

where we have denoted by r+ the radial position of the horizon. In the r — 
R plane, the metric (|42|l becomes then regular on the horizon i? = if the 
complexified time r has a period /?; see e.g. p31 ESI IMI 15^ . 

Remark 3.1 In spacetimes with multiple horizons Hi there is no global notion 
of time-periodicity. Basically, the neighborhood of each horizon can be covered 
by a system of coordinates in which the metric is regular on the horizon and 

*A wider class of metrics will be considered in the next section. 
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which provides us with a well-defined horizon temperature Ti = l/Pi- However, 
in the spacetime regions covered by more than one coordinate system, different 
notions of temperature exist which in general are incompatible among them; see 
j59|. Despite of this, the definition (|41|l of temperature still has a local meaning 
on each single horizon. ■ 

Remark 3.2 Let us now consider a surface S of constant time with an outer 
boundary denoted by B. Let us also suppose that S intersects a certain number 
of horizons Hi with intersection surfaces denoted by Hi. The (variation of the) 
internal energy contained in St is given by 6xE{£^, B). Since we are dealing 
with stationary spacetimes and since the surface B is homological to the union 
of all Hi (i.e. B Ui Hi = is a homologic boundary) expression (|24|l holds 
then true. Therefore the variation of the internal energy receives a contribution 
from each cross section Hi: 

5xE{i,B)^Y.^xE{^,Hi) (44) 

i 

Each Hi being endowed with its own energy content and with its own temper- 
ature Ti we expect it features also an entropy contribution Si according to the 
local first law 

5xE{(„H,)=T,5xS^ (45) 

Basically, for a spacetime with multiple horizons, where a global notion of tem- 
perature does not exist (namely, the spacetime does not describe a thermody- 
namical system in equilibrium) , there cannot hold a global first law. Neverthe- 
less many local first laws hold according to (|45|l . 

We stress that expression (^51 is calculated directly on the horizon indepen- 
dently of what exists outside it. This is a rather desirable property. Indeed even 
if the whole gravitational system is not in thermal equilibrium it is nevertheless 
expected that an observer located near the horizon should be able to measure 
physical properties of the horizon itself without knowledge of spacetime regions 
far away from him. 

We also point out that the expression as it stands, is suitable for appli- 
cations to a broader class of horizons, e.g. isolated horizons 0, the geometric 
characterization of which is intrinsically independent on the geometric proper- 
ties of the spacetime surrounding them. In the isolated horizons framework, the 
whole spacetime is not even required to be stationary. Rather, it is just assumed 
that no flow of matter/energy crosses the horizon. Under these conditions the 
horizon is isolated and it can be considered as a thermodynamical system in 
equilibrium. Therefore the laws of thermodynamics can be generalized also to 
this class of horizons; see jHlini- In it was indeed shown how the first law can 
be geometrically reproduced by means of H45|l . 

We finally point out that if we had made use of the definition (|33f) of internal 
energy the above construction would have failed to work since £{g,Hi) — 
(being both the lapse and the shift vanishing on Hi). 
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In the purely Lorentzian sector, we recover the Lorentzian counterpart of 
expression (|42l) by considering the class of metrics [5^1 which are: 

i) static in the given reference frame, 

ii) have vanishing lapse on some (compact) 2-surface Hi defined by iV^ = 
(surfaces corresponding to removable singularities of spacetime) and 

iii) non vanishing derivative 'VaN'^ 7^ on Hi. 

Near each horizon Hi we can approximate the lapse function with its Taylor 
series truncated to the first order, obtaining in this way the metric: 

32 /'27ry , . 



g = -R^ dt' + dR' + r\R)dn (46) 

where T = 1//3 is given by (|41|l while r{R) corresponds to H43|l . Namely: the 
zeroes of N'^ determine the horizon radii r_|_ while the radial derivatives of N'^ 
fix the temperatures /?. When we consider variations Sxg along a generic vector 
field X, both parameters r+ and /3 are varied. 

For the class of metrics H46f) the boundary condition (|35|l calculated on the 
horizon {R = 0} becomes: 

2 

TTT" 

^^^5IinH = -^^P = ^r\5T (47) 

Therefore, the (weak) Neumann boundary condition ^^u^^^'^\h = is sat- 
isfied if <5/?|^^o = 0. Namely, the Neumann ensemble is locally realized by the 
set of metrics admitting a causal horizon inside the surface S with a fixed tem- 
perature T = i. Hence, for vector fields X which describe variation along the 
Neumann ensemble, from H45|l we obtain: 

5^M](TS) = T5^M]{S) = 6^N]E{^, B) = 5^m] (^J^ C/Kom(o) (48) 

where in the first equality we made use of the property d[i\f^{T) = while the 
last equality comes from expression (|36|) . The above relation can be formally 
integrated as: 

TS= f UKon.{0+Q0 (49) 

where Qq is a constant of integration and corresponds to the Komar potential 
calculated for a (background) solution gg lying in the integral curve of the vector 
X. It corresponds to the ground state for the calculation of the gravitational 
heat TS. From now on we shall set the constant Qo equal to zero so that the 
reference solution corresponds to the one for which the horizon H shrinks to 
zero, i.e. the solution for which — (since our analysis is only local we do 
not forbid the background solution to have other horizons elsewhere). 
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It is now easy to verify that the entropy, as defined by (|49ll . corresponds to 
one-quarter of the horizon area Ah- Indeed, when we approach the horizon Ti. 
the metric assumes its form (I46|l and expression (|37f) becomes: 




provided that the orientation of the unit normal is chosen so that iVn^a^/27r 
turns always to be positive.^ 

We stress that the definition of entropy as it arises from (|49|l and l|5Q|l 
is clearly observer-dependent. Indeed, conditions i)-iii) listed above may be 
rephrased as follows. Given a set of observers evolving with 4- velocity , let us 
choose an observer-adapted system of coordinates in which the observers them- 
selves are at rest. In this reference frame the metric coordinate singularities 
correspond to homological boundaries and, according to each boundary 

gives rise to its own contribution to the total entropy. However, from a physical 
point of view, the metric coordinate singularities correspond to regions of inac- 
cessibility for the observers (one-way membranes). Entropy turns then out to 
be associated to the observer unseen degrees of freedom, namely, to the regions 
which the observers have no physical access to. According to [^^1 we might 
be tempted to state that entropy measures the information content beyond the 
hidden regions and that such an information is related to one-quarter of the 
area of the surface enclosing the observer-unaccessible regions. 

Let us now end this section with some further thermodynamically-inspired 
consideration. Let us consider again expression (|25|l and let us evaluate it on 
the cross section Hi of a causal horizon. According to H49|l it can be rewritten 
as: 

SxE{^, H,) = 5x{TrS,) + f ^g^"^ (5u[f, C^ds^f^ (51) 

From a thermodynamic point of view, having formally identified E{^, Hi) with 
the internal energy of the gravitational system and the Noether charge with the 
TS contribution, it is tantamount to formally identify the second term in the 
right hand side of (|51|l with the variation of the gravitational free (or Helmholtz) 
energy F{^,Hi): 

5xF{i,Hi) := / :/lg^^- 5u%i''Us^p = - f -i^JWcPx 

= I UcaM^.X) (52) 

Accordingly, expression (|51|l turns out to be nothing but the Gibbs-Duhem 
formula ^]^: 

SxE{^,H,) = 5xiT,S,) + SxFit Hi) (53) 

^For non compact horizons, e.g. Rindler horizons ^9 , formula I5UI gives a finite value 
only if the domain of integration extends up to a finite spatial region. Despite of this, the 
entropy for unit of area is always a finite constant and it equals 1/4. 
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Trying to further enhance the formal analogy between gravitation and thermo- 
dynamics we observe that, in the definition (|45l) of the first law, we are taking 
into account variations for which no "mechanical work" is done. Namely, all the 
variation of energy 6E among nearby solutions corresponds to "heat" exchange 
T5S. For such (reversible) tranformations, thermodynamics predicts that the 
free energy F must obey the law ^ = — S'. This is exactly in agreement with 
the identification (|52() and with equation (|47|l . For instance, if we consider two 
nearby Schwarzschild solutions (of mass, respectively, M and M + 6M) the 
variation of energy SE turns out to be SE = SM (see appendix ^ and it is 
equipartitioned into variation of "heat" and variation of "free energy" according 
to SM = S{M/2) + S{M/2). This splitting corresponds, from a thermodynam- 
ical point of view, to the splitting of SE = TSS into TSS = 6{TS) - SST. 
The free energy contribution S{M/2) therefore corresponds to the —S6T part 
of the splitting. It does not correspond to mechanical work (no energy can be 
extracted from the hole in any classical physical process; see e.g. |SZ1) but it 
is due to the fact that the transformation is not isothermal (if we change the 
mass of a Schwarzschild black hole we also change its temperature according to 
T = I/SttM). Indeed -SST = -AttM^ 5{l/8nM) = 5{M/2). 

Further hints corroborating the viability of our definitions will be exhibited 
in the appendix El where the analogies of the formalism developed so far with 
existing literature on the subjet will be analysed. 

4 Generalized Entropy 

Following the guidelines of [SJj and adapting Mann's ideas to our formalism, 
we shall now try to generalize the identification of the Noether charge with 
gravitational heat to a wider class of situations. We nevertheless stress that our 
approach differs considerably from that of Even though the latter approach 
turns out to be more physically motivated we believe that our definitions are 
better suited to handle specific apphcations, e.g. the AdS-Taub Bolt solution, 
where comparable results can be obtained with considerable less efforts. 

Let us then consider a stationary spacetime M . Let ^ = 9t be the Killing 
vector field and let us then denote by UiTii the union of the fixed point sets of ^. 
Any T-Li is an obstruction to foliating spacetime with surfaces of constant time. 
Therefore we can consider a surface S of constant time the boundary of which is 
formed by an outer boundary B together with the union of all the cross sections 
Hi of Tii. If a temperature Ti can be defined for any 7?^, following step by 
step the same arguments outlined in Remark l3. 21 we expect that each Hi gives 
rise to a contribution TiSSi to the internal energy and therefore a contribution 
Si to the total entropy. Notice however that now our definition encompasses a 
broader class of solutions, namely the ones with any kind of fixed set points of 
the Killing vector. For Lorentzian spacetimes, it includes any compact causal 
horizon but also non compact horizons, such as Rindler horizons or Lorentzian 
CCBH spacetimes (constant curvature black holes |22j). When dealing with 
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Euclidean gravity, bolts, nuts as well as Misner strings (HHI are allowed to exist. 

To calculate the contributions to the entropy arising from each singularity 
Hi, it is a common procedure to start from the (local) relation: 

TA = E{Hi) - F{H,) (54) 

relating temperature and entropy with internal and free energies of each singu- 
larity Hi. The idea to make use of the above relation dates back to Gibbons 
and Hawking |37j . They identified the internal energy with the Hamiltonian 
and the free energy with the gravitational action, as defined in formula (|1(J3|I of 
appendix In general both the Hamiltonian and the action, when calculated 
on a given solution, give an infinite result. Therefore one has to introduce a 
background solution, suitable matched with the original solution, and subtract 
off from H54(l the background action and Hamiltonian. In this way one obtains 
a finite result which provides us the difference of entropy between the solution 
and its background. 

Another strategy which avoids the problems related to any background fixing 
procedure has been developed in |54| . It consists in modifying the gravitational 
action, i.e. the gravitational free energy, by suitably adding boundary counter 
terms coming from the conjectured AdS/CFT correspondence. Gravitational 
internal energy is then obtained, via a Brown- York procedure, from the modified 
action functional. 

We shall instead face up to the problem following the technique of section 
13 First of all let us notice that if we consider the infinitesimal version of (|54|l : 

5{T,S,)^5E{H,)-5F{Hi) (55) 

the problems of fixing a suitable background is postponed, in specific appli- 
cations, up to the end of calculations where the background will resort as a 
constant of integration. This implies a higher practicality since it is not nec- 
essary to select from the very beginning the reference solution and to carry on 
calculations with it. Indeed, according to our definitions H38|) and (|52|) . the 
relation H55|l becomes: 

,5(T,50 ^(^ly^ C/Kom| (56) 
which, after integration, leads to the result 

T,S, = / C/Kom + Qo (57) 

Now, it is the constant Qo which encodes the information relative to the back- 
ground. Different choices of Qo allow to reproduce many results exhibited else- 
where. For instance, Qo can be set equal to zero for compact horizons. This 
choice, see (|50|) . leads indeed to the one-quarter area law for the entropy. Other 
choices are nevertheless possible and strictly necessary for non compact hori- 
zons, where a background subtraction term is necessary to obtain a finite result. 
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To this end, we recall that, given a variation 5x induced by a vector X in the 
space of solutions, the background has to belong to the integral curve of X 
through the original solution. This means that the given solution has to be, in 
the suitable sense, continuously deformable into its background, i.e. the solution 
and its background have to lie, at least, in a path-connected region in the space 
of solutions. This consideration suggests how the background has to be chosen 
in many of the applications. The example below will help to clarify this point. 



5 AdS Taub-Bolt solution 

Let us consider the class of metrics of the kind |18l I54| : 

g = V(r) (dr + 2Ncosed(l)y + —— + (r^ - N^)(d0^ + sin^ 0d(l)^) (58) 

V{r) 

+ 7V2 _ 2mr + (r* - 6N^r^ - 3N^) 
V[r) = — TO, N constant 

Independently on the parameter to, they are all asymptotically locally AdS |43) : 

g = Voc.(r) (dr + 2N cos ed(t)f + + (r^ - N^){de^ + sin^ edcf)'^) 

V^{r) = ^+[l-—j+0{r-') (59) 

The asymptotic boundary is a squashed 5*^, rather than x |43[ I56| . and 
the constant N parametrizes the squashing. When we apply the definition H22|l 
to the metric (|58|l . variations 5 are meant inside the class of solutions with 
asymptotic behavior H59f) . After some calculations we obtain: 

SxEi^,SR)^j^J^Sm (60) 

where Sn is a surface of constant radius R. In the limit i? ^ oo we have 

SxEi^,S^.)^dm (61) 

so that we can refer to the parameter m as the total Energy/mass of the solution. 
Notice that, even though the quantity SxE is known to be an homological 
invariant (see equation (|24|) ') formula (|60|l depends on the radius. The reason is 
easily understood. 

One of the fixed points set of the U{1) isometry r is the surface r = vq where 
To is a solution of V{rQ) = 0. The equation V{rQ) = allows to determine the 
"mass" parameter to as a function of tq. From 1)58(1 we indeed obtain: 

TO(.o) = r^o^jP-m^H^NHP-m^) ^^^^ 

Zt Tq 
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The fixed set points r = tq is called a bolt and its area is equal to 



A{ro) = SttN 



(63) 



Moreover, the metric (|58() has another coordinate singularity, i.e. another break 
down of the foliation in constant time surfaces, due to the presence of a Misner 
string. It is a two dimensional coordinate singularity running along the z-axis 
(i.e. 9 = 0,tt) from the bolt out to infinity (see ISEl)- The Misner string was 
recognized in |56| to be a removable singularity provided we fix the period (3 
of the Euclidean time to be = SttN. In addition, the absence of a conical 
singularity at the bolt r — rQ determines the constraint |43|: 



f3 = SttTV : 



47r 



V'iro) 



(64) 



The above equation admits only two solutions, i.e. rg = and tq = r„, where: 



12N 



= N 



(65) 
(66) 



for the values N < '''^^j^^-'' corresponding to a real square root. This two 
solutions are known, respectively, as the Taub-bolt and the Taub~NUT metrics. 
By substituting tq = rn into l|63|l . we easily verify that the NUT has zero area. 

From the above discussion it follows that a surface enclosing the bolt (or the 
NUT), i.e. where tq = r^, (or, respectively, tq — rn) is not homologous 

to spatial infinity. Owing to the presence of the Misner string we also have to 
consider two cones Ci ~ {9 = e,ro + e < r < 00} and C2 = {9 = tt — e,ro + e < 
r < 00} wrapping around the z-axis from the bolt r = tq up to spatial infinity. 
Spatial infinity is then homologous to the union S = S^ra+e) U Ci U C2 . From the 
homological property (|24|l we known that each term in S gives a contribution 
to the total energy H61|l and, accordingly, to the total entropy. 

We can calculate each separate contribution to the entropy following the 
definition (|57|l . We obtain 



TS. 



for the bolt contribution and 



lim 

£^0 



T 



(67) 



r=ra+<i 



{C1UC2} 



lim 



t/Kom + / UKom 
=e Je=Tr-e 
2 /\r2/q)\r2^ I ^3 ;2^ , ^("^ \ 72\ 



(68) 



= lim <, -rr^r H ttttttt^ 2\ ^ 



for the Misner string contribution (the details of the calculations follows the 
ones described in |24[ I34| for the asymptotically locally flat case P 00). The 
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total entropy Srg is given by the sum of and 



Sro = — i lim 

4 r— »oo 



+P ■ Qo (69) 

where Qo is the constant of integration; see H49|l . If we do not take into consid- 
eration the constant of integration, the total entropy diverges due to the 
infinite contribution coming from the Misner string (while it reproduces a finite 
value in the limit P — > oo). Notice however that the divergent part turns out 
to be independent on tq. Therefore if we consider the entropy of the Taub-bolt 
solution relative to the Taub-NUT metric we obtain a finite result and, in the 
meanwhile, we get rid of the constant of integration which is common for both 
the solutions. Indeed both Taub-bolt and Taub-NUT can be joined, in the 
space of solutions, by a curve 7 the points of which all share the asymptotic 
form H59|) . This can be done by just varying rg from r„ to ri,. The constant 
Qo corresponds to the Komar integral calculated on an (arbitrarily) fixed point 
go on 7 ^. The relative entropy of Taub-bolt with respect to the Taub-NUT is 
then given by: 

AS = Sr,-Sr = ^ \H + rliP - 12iV2) + 2nN{6N^ - P) + N\l^ - ZN^)] 

(70) 

for both the values = . We remark that the entropy H70|l satisfies the first 
law 

T6{AS) = Sm{n) - (5m(r„) (71) 
and it also agrees with the results of [TS|. 



6 Conclusions 

We finally point out that the thermodynamic analysis we have performed en- 
tirely develops from the definition H22|l of variation of energy. We remark that 
this definition, which perfectly agrees with many covariant definitions obtained 
elsewhere (see, e.g. |EI UHl El ESI E])) was calculated directly from the 
Hamiltonian equations of motion and therefore it does not depend on the spe- 
cific gravitational Lagrangian L inside the variational cohomology class [L] we 
choose to represent the system, where elements of [L] differ each other only 
for the addition of divergence terms. About this, let us again point out that, 
even if throughout the paper we have dealt with the Noether charge, nowhere 
in the text we made use of Noether theorem, which is instead sensitive to the 
representative chosen inside [L]. The splitting of SE into 6{TS) and 6F depends 
therefore just on the field content of the theory (i.e., Euler-Lagrange equations) 

^In particular, diverging terms are eliminated from 1169 1 for any choice of the background go 
in 7 (and none of them corresponds to Qo = 0)- Nevertheless background choices other than 
the Taub-NUT solution turn out to be, as far as we know, meaningless and rather artificial. 
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and not on a given Lagrangian description inside [L] . Only if we try to formally 
integrate these functions of state with the appropriate boundary conditions we 
restore the Lagrangian framework and it can be shown (see appendix^ that 
standard (and widely accepted) results are reproduced in this way. Nevertheless 
we stress again that in specific applications it is more convenient for calcula- 
tions to consider just variations of physical observables without performing any 
a priori integration. In this way complicated calculations related to background 
choices and background matching conditions are avoided in a first approach. 
The background will instead come into play just at the end, in a more manage- 
able fashion, as a "constant of integration" . 

As future developments, we believe that the ideas suggested in this paper 
are worth being extended to deSitter and Anti deSitter Einstein-Gauss-Bonnet 
gravity in order to study the possibility of negative entropy; see Moreover, 
the analogy between Noether charge and gravitational heat deserves to be ana- 
lyzed in processes of collapse of a star into a black hole in order to understand 
if a sort of gravitational heat can be defined prior to collapse. 

We would like to thank S. D. Odintsov and R. B. Mann for having addressed 
our attention on these future fields of investigation. 
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A Spherically Symmetric Solutions 

Let us consider a static spherically symmetric solution of the kind: 

dr'^ 

g = -f{r)dt^ + — +rHn (72) 

fin 

where dVl — dO"^ + sni{d)^ dcfP' . The explicit form of f{r) is dictated by Einstein's 
field equations and we shall analyze it later on. 

The variation of the energy (|22|l contained inside a region St of constant 
time bounded by a a two-sphere Sr, i.e. Bt := {r = i?}, turns out to be: 

5xE{^,Sr)^ ( d^xN5{^e) = -^8f\^ (73) 



21 



where £^ — dt- The Noether charge (I37|l . apart for the constant of integration 
which is set equal to zero throughout this section, becomes: 



Q{^,Sn)^^\nR)\ 



(74) 



having denoted with a prime the derivative with respect to the r-coordinate. 
The absolute value arises because when /' > 0, e.g. on black hole horizons, 
the normal to the sphere is chosen to point outward, along the direction of 
increasing r. When /' < 0, e.g. on cosmological horizons (see below), the unit 
normal is instead chosen to point inward toward the direction of decreasing r. 
In this way the term " Jj''^ = in the Noether charge (|T7|l turns always 

to be positive. 

Dirichlet boundary conditions for (|73|l have been extensively analyzed in 
[Tl] . We instead point out here that (weak) Neumann boundary condition 
corresponds to: 



/ ^^JW^ ^~^{R5r{R) + 25f{R)} = 

^ Sr 



(75) 



It is therefore easy to check that Neumann variations (5[Ar](9(^, Sr) of the Noether 
charge (|74f) give rise to the variation of energy H73|l . No anomalous factor arises 
in this case! 

Now, let us suppose that /(rj — for some r.;, (i — 1,2,...) and that 
f'{ri) 7^ 0. We define the temperature T,; of each horizon r — Vi according to 
gll), i.e.: 



On the horizon, (|75|l is equivalent to 6Ti — and the Neumann energy H74|l 
becomes: 



fin) 



4tt 



A, 



(77) 



where Ai is the area of the horizon's cross section. Therefore the above expres- 
sion has the TiSi form if the entropy is identified with Si := We stress that 
so far we have not required the solution to be a vacuum solution (possibly with a 
cosmological constant). Therefore the result (|77|l holds true for any spherically 
symmetric static solution. 

Further considerations require to consider on separate grounds the cases in 
which: 

i) the metric (I72|l is a solution of Einstein's equations with a (positive) cosmo- 
logical constant; 

ii) the metric is a solution of Einstein's equations with external matter fields. 



i) Schwarzschild— deSitter solutions 

If we set 



2to 

f{r) = 1 — 

r 



(78) 
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in the metric (|72|l we obtain the Schwarzschild-deSitter metric. The cosmolog- 
ical constant is A = 3//^ and m is the mass parameter. If m = we recover 
deSitter spacetime while Schwarzschild solution is obtained for A 0. 

Provided that m < I/^/TI the solution has two positive roots, r+ and r++ 
(r+_|_ > r+) corresponding, respectively, to the black hole horizon and the cos- 
mological horizon. From (|78|l we have: 

for both the values ro — r+ and = 

When r ranges from r+ up to r-|-+ the vector field ^ = dt is timelike and 
it can describe the velocities of observers located between the two horizons (of 
course this is just one possible choice; see for another normalization of ^). 

Expression H73|) becomes: 

SxEi^, Sr) = Sm (80) 

and does not depend on the radius of integration, as expected from H24|l . 
The temperatures on and r++ are, respectively: 

f'{r+) 
An At: \r+ P 

where the minus sign in H82|) arises since the surface gravity ^^2^^ negative 
on the cosmological horizon. According to (|77|l the integral of the Komar charge 
allows to obtain the one-quarter law for both the horizons ro =7"+, ?"++: 

5*0 = nrl (83) 

Notice from H79(l . H81(l and H83(l that, on the black hole horizon, the relation 
5m = T+(55+ holds true. Owing to H8U|) this relation corresponds to the first 
law 

5E{£„ r+) ^6m^ T+SS+ (84) 

On the cosmological horizon, on the contrary, it holds true the relation 5m = 
—T++5S++ so that we expect (5iJ(^,r++) to be equal to —5m if we want to 
recover the first law for the cosmological horizon, too. This result is only appar- 
ently in contrast with H80|l . Indeed, in obtaining H80|l the spacelike unit normal 
was assumed to point in the direction of increasing r. Nevertheless on the 
cosmological horizon we have to change the sign, i.e. 

5xE{S,,r++)^-5m (85) 

if we want the unit normal to point toward the observer. In this way the 
expected value is recovered and the first law is satisfied. Notice that a nega- 
tive value for the variation of energy contained inside a cosmological horizon is 



23 



physically expected; see ^1 |H21 ■ Indeed, if a test particle of mass Sm is thrown 
across the cosmological horizon from the interior region, the area of the horizon 
increases but we expect that the energy enclosed by the horizon should decrease 
according to (|85|l . 

Notice that, even though the absolute value \dE\ — Sm between two nearby 
Schwarzschild-deSitter solutions is the same when computed on both the hori- 
zons, the energies of the horizons i?(^,r++) and i?(^,r+) can be nevertheless 
different from each other (even though, in absolute value, they change for the 
same amount during the variation). From (|84|l we obtain: 

= Jo ^ ^ "^'^^ ^ 

where the normalization has been reasonably fixed so that the energy vanishes 
when the horizon mass vanishes, i.e. the background spacetime is the deSitter 
solution. Notice that, for — > oo the usual mass ^ of Schwarzschild black hole 
is recovered. 

According to H85|l , for the cosmological horizon we obtain: 

E{C,r++) = -mir++) + Eo (87) 

In this case different choices are available for the constant of integration. For 
instance, if we set i?o = we fix again the deSitter background. Another pos- 
sible choice is Eq — -^=- In this case E{£_,r++) E{£^,r+) in the Nariai limit 

m I/V27 when r+ and approach to each other. We refer the reader to 
[T!?] for some deeper insight into the subject (see also I85j). 

ii) External matter fields 

If we set 

/(r) = 1 - !± - 1 r e{f)f^df (88) 

in the metric (|72(l we obtain a solution of Einstein's field equations with external 
matter characterized by a density 

Pir) = (89) 

radial pressure 

p{r) = p{r) (90) 

together with angular pressures Tg = = -j^, the explicit form of which is 
of no interest here; see 1301 EH|- The function (|88|l has been normalized in 
order to have an horizon at r+ with temperature: 
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We however do not exclude that /(r) can be zero for other values of r. 
Expression H73() now becomes: 



R 

1 



e{f)f^df } (92) 



for a generic sphere of radius R and 

6xEi^,r+)^s{^-±} (93) 

for the horizon. 

If we consider variations which keep the radius R fixed in ()92fl . we easily 
obtain: 

SxEi^, Sr) = SxEi^, r+) - ie(r+)r^<5r+ (94) 

Since from (|91|l we have: 

6xE{^, Sb) = T+SS, S - nrl (95) 

by comparing (|94|l and H95|) we obtain the first law: 

SxE{^,r+) = T+SS+pJV^ (96) 

where p+ is the radial pressure H90() on the horizon, while SV^ = 47rr^^r+. 
Hence p_^_SV^ is the work term given by the product p^Anr"^ of the radial force 
multiplied by the radial displacement Sr^. 

Therefore expression H96|l constitutes the generalization of formula H45|l in 
presence of external matter. Namely, the amount of energy to virtually increase 
the radius of the horizon of a quantity Sr^ is partitioned into a part T^6S which 
comes from the increase of entropy (due to the expansion of the horizon area) 
and a work term p^5V^ which has to be spent to contrast the external pressure, 
see jSHl- 



B Comparison with the semi— classical approach 



Let us now analyse the interplay between the formalism developed so far and the 
existing literature on the matter. First of all we stress that the identification of 
the Komar energy H37(l with the gravitational heat TS ensues from the heuristic 
definition of entropy given in (|45|l . which a priori has no direct physical ground. 
The identification is nevertheless in agreement with semi-classical statistical 
approaches, based on path integral techniques, where entropy is identified by 
the value of the microcanonical action functional (evaluated on a -complexified- 
stationary black hole solution; see [111 1141 US) ). 



25 



Indeed, let us consider a closed surface B surrounding the horizon H such 
that BUH is the homologic boundary dT, of a region S. We can rewrite equation 
as^: 

SxiTS) = 6xE{^,H)-SxF{i,H) 

121 



SxE{^,B)-6xF{^,H) 

t/Kom| + C/cadm(C, ^) 



B 



Sx< / t/Kom+ / ii\L\ (97) 



where L = ^ '^•^ Hilbert Lagrangian. In the last equality we 

made use of Stokes' theorem together with the following relations which hold 
true on-shell for stationary solutions: 



^ dC/cADM(e,^) (98) 
From l|97(l we obtain (apart from a constant of integration): 

S - /3|^C/Kom+ / «d^l 



J dtj C/Kom+ J dtj 



B J JT. 

dt/\UKom+ I L = I.,n (99) 

where Im is the microcanonical action functional for a region D — P x T, oi 
spacetime with outer boundary B — (3 x B ; see, e.g. ^103]. Hence definition 
H49() together with the relation (|53|l exactly corresponds to Brown- York's defi- 
nition of entropy when stationary black holes are considered. 

We also stress that the identification S = /3 UKom exactly agrees with the 
definition of entropy |^ : 

S^-[ ^d'^xMn'^a^, (100) 



^For notational convenience, from now on, we shall suppress the index i. 
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which Padhmanabhan has given for static spacetime, as it is can be easily in- 
ferred from (|37|l . Nevertheless the definition (|52|l of free energy differs from 
Padhmanabhan's definition j59j in which free energy is identified with the Hilbert 
Lagrangian. Our definition is instead more tied with the Trace-K action func- 
tional By making use of the relations H98|) . definition (|52f) can be indeed 
conveniently rewritten as: 

-SxF{^,H)= f t^lSL- f Ucaum{^,X) (101) 

JS J B 

If we consider the Dirichlet boundciry condition Sx^fii^ L = on the outer 
boundary B (which in many applications can be identified with spatial infinity), 
namely, if we consider the class of metrics which approach a fixed background 
go on B, the above expression can be integrated. From H28() we hence obtain: 

~F{(,H)= [ i^L-- [ y/^d^xN{e-Qo) (102) 

Jn ^ J B 

After multiplying the above expression by dt and integrating over the time 
interval [3 we finally obtain the trace-K action functional Ig-. 

Ig:= [ L-- [ ^d^xN (6 - Go) = -/3 • F(C, H) (103) 
Jd Jb 

(which coincides, in the hypotheses we are dealing with, with the first-order 
covariant Lagrangian for General Relativity; see [211 )■ We recall that, in a 
semi-classical path integral approach, Ig corresponds, in the Euclidean sector, 
to the logarithm of the thermodynamic partition function; see jl4l 
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